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ABSTRACT: We investigate non-thermal gravitino production after tribrid inflation 
in supergravity, which is a variant of supersymmetric hybrid inflation where three 
fields are involved in the inflationary model and where the inflaton held resides in 
the matter sector of the theory. In contrast to conventional supersymmetric hybrid 
inflation, where non-thermal gravitino production imposes severe constraints on the 
inhationary model, we find that the “non-thermal gravitino problem” is generically 
absent in models of tribrid inflation, mainly due to two effects: (i) With the inhaton 
in tribrid inflation (after inflation) being lighter than the waterfall held, the latter 
has a second decay channel with a much larger rate than for the decay into gravitinos. 
This reduces the branching ratio for the decay of the waterfall held into gravitinos. 
(ii) The inhaton generically decays later than the waterfall held, and does not produce 
gravitinos when it decays. This leads to a dilution of the gravitino population from 
the decays of the waterfall held. The combination of both effects generically leads 
to a strongly reduced gravitino production in tribrid inflation. 
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1 Introduction 

Cosmic inflation is a successful paradigm for the early universe. Among the various 
types of inflationary theories proposed, “hybrid-like” models [1] offer the attractive 
possibility to connect them to a particle physics phase transition around energy scales 
Mqut of Grand Unification. In such models, inflation proceeds along a nearly flat 
valley in held space with large vacuum energy. When a “critical value” is passed, 
inflation is ended by a second order phase transition where the so-called waterfall 
held moves quickly towards its minimum. 

In the “conventional” realisation of hybrid inflation in supersymmetry (SUSY) 
or supergravity (SUGRA) [3-5], the waterfall held can be associated with the spon¬ 
taneous breaking of some symmetry around Mq\j t, however the inflaton held is 
necessarily a singlet under any symmetry (apart from possible R-symmetries). It is 
therefore somewhat decoupled from the rest of the particle theory. 

This can be overcome in tribrid inflation [6, 7]. In tribrid inflation, three fields 
are involved in the inflationary model: the inflaton held, the waterfall held, and a 
so-called “driving held” (which however does not contribute to the dynamics during 
and after inflation). In contrast to “conventional” hybrid inflation, the inflaton held 
can now reside in the matter sector of the theory. One can, for example, identify 
it with one of the right-handed sneutrinos [8] or with a D-flat direction of fields 
in representations of the gauge symmetry [9]. This allows for particularly close 
connections between inflation and the underlying particle physics theory (see e.g. 
[ 10 , 11 ]). 

Inflationary models of the early universe face various constraints: For instance, 
they have to be consistent with the data on primordial fluctuations extracted from 
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CMB observations, which requires a comparatively flat potential during inflation. 
After inflation, the models have to successfully reheat the universe and produce a 
thermal bath of Standard Model particles. In addition, not too many long lived 
cosmic relics should be produced. One example of the latter problem in the context 
of supergravity is the so-called “gravitino problem” [12]. 

Gravitinos are the superpartners of the gravitons in SUGRA. Since they interact 
quite weakly (with strength of gravitational interactions) they are typically rather 
long lived and can be problematic in various respects: In particular, if they are stable, 
they can easily overdose the universe. When they are unstable and decay after BBN 
(or around the time of BBN), their energetic decay products can destroy the produced 
light elements. But also if they decay before BBN, which is typically the case when 
they are heavier than about 30 TeV, they impose significant constraints. The decay 
of each gravitino (assuming R-parity conservation) in general produces one stable 
lightest supersymmetric particle (LSP), which would count towards the dark matter 
(DM) relic abundance. The constraint on the number density of gravitinos in this 
case comes from the requirement not to overproduce DM. 

Gravitinos can be produced thermally, in the thermal bath after reheating [13- 
15], and non-thermally from the decays of the fields involved in the dynamics after 
inflation [16, 17]. The thermal production imposes constraints on the reheat temper¬ 
ature. For instance, assuming a gravitino with mass above 30 TeV (with conserved 
R-parity) and a LSP mass of about 100 GeV, this bound is roughly T R < 0( 10 10 GeV) 
(see, e.g., [18]). For much smaller T R , the thermal production can be neglected in 
this example case. On the other hand, it is known that for the case of “conventional” 
SUSY hybrid inflation the non-thermal production of gravitinos provides a potential 
problem [17]. 

In this work, we therefore investigate non-thermal gravitino production after 
tribrid inflation and compare the results to the case of “conventional” SUSY hybrid 
inflation. We find that the “non-thermal gravitino problem” is generically absent 
in tribrid inflation, mainly due to two effects: (i) With the inflaton being lighter 
than the waterfall field, the latter automatically has a second decay channel with 
a much larger rate than for the decay into gravitinos. (ii) The inflaton generically 
decays later than the waterfall field, and does not produce gravitinos when it decays. 
This leads to a dilution of the gravitino population from waterfall field decays. The 
combination of both these effects generically leads to a strongly reduced gravitino 
production in tribrid inflation. 

In the next section, Sec. 2, we will review “conventional” SUSY hybrid inflation 
and tribrid inflation, with a focus on the post-inflationary dynamics. In Sec. 3, we 
discuss the non-thermal gravitino production from a generic scalar field at the end of 
inflation, followed by the calculation of the gravitino yield in “conventional” SUSY 
hybrid inflation in Sec. 4. In Sec. 5, we calculate the gravitino production in tribrid 
inflation and discuss our results. Our conclusions are presented in Sec. 6. 
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2 TTibrid inflation and conventional SUSY hybrid inflation 


There are two classes of supersymmetric hybrid-like F-term inflation models: “Con¬ 
ventional” F-term SUSY hybrid inflation [3] and tribrid inflation [7]. Both have in 
common that there is an inflaton field responsible for the dynamics during inflation 
and a waterfall field responsible for ending inflation. The end of inflation may be 
associated with a particle physics phase transition, which opens up the possibility to 
embed inflation in particle physics models such that the combined theory becomes 
more predictive and testable. 

“Conventional” F-term hybrid inflation: When SUSY (or SUGRA) hybrid 
inflation is discussed, the model is usually described by the superpotential 

W = K<t){H 2 - M 2 ), (2.1) 

combined with a (at least nearly) canonical Kahler potential [3-5] . The superfield 0 
contains the inflaton field as a scalar component. For simplicity, we have assumed 
H to be a singlet, but H 2 can be replaced by HH where H and H belong to the 
conjugate representation of a gauge group. We will use the same symbols for the 
superfield and the scalar component throughout the paper. 

In conventional F-term hybrid inflation, 0 is a singlet*, 0 is the so-called “driving 
field” which generates the potential for H by its F-term \F ( p\ 2 , but it is usually 
somewhat disconnected from the rest of the particle physics theory. During inflation, 
(H) = 0 and 0^0 (false vacuum), driving inflation. The tree-level inflationary 
potential Vo = ft 2 M 4 is corrected by quantum corrections (and Kahler potential 
corrections) generating the desired slope of the potential. Without loss of generality 
we consider real k and M. With k ~ 0.1, an appropriate matching of the amplitude 
A s of the scalar fluctuations requires M ~ 5-10 15 GeV. On the other hand, an exactly 
canonical Kahler potential would be inconsistent with the PLANCK observation of 
the spectral index [2]. This tension can be resolved, for instance, with higher order 
terms in the Kahler potential for the inflation field [19], or by considering the effects 
of soft supersymmetry breaking terms [20]. 

H is the waterfall field which ends inflation by a second order phase transition 
when 0 rolls below a critical value where the mass of H gets tachyonic. After this 
“waterfall” (assuming homogenous fields) both fields, the inflaton and the waterfall 
field, oscillate around their minima. When the SUSY breaking effects are neglected, 
the minima are at (0) = 0, and (H) = M. Including SUSY breaking effects after 
inflation, the waterfall field and the inflaton field mix almost maximally and are 
almost degenerate in mass, 

0 ± =-b(*±ff), (2.2) 

with m 2 ± ~ 4 k 2 M 2 . Let us now turn to tribrid inflation, which in several respects 
has quite different features. 

*It may only be charged under some R symmetry. 
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“Tribrid” inflation: In tribrid inflation [6, 7], the inflaton can be associated 
with a matter field of the underlying particle physics theory. Here, for simplicity, we 
consider the case that it is a singlet as well (e.g. one of the right-handed sneutrinos 
in extensions of the MSSM to explain the observed neutrino masses) but it may also 
be a D-flat direction of gauge non-singlet matter fields as discussed in [9, 21]. There 
are several variants of tribrid inflation in supergravity [6, 22, 23], but we will focus 
here on a simple representative where (in Planck units) the superpotential has the 
form! 

W = k X {H 2 - M 2 ) + A H 2 cj) 2 . (2.3) 

The driving field is now called A". It is not the inflaton and only provides the 
potential for the waterfall field H by ist F-term F x | 2 • By a term 7 _\-|A"| 4 in the 
Kaliler potential K with yw < — 1 /3 it will have a mass larger than the Hubble scale 
during inflation and thus stays at ( X) = 0, not affecting the inflationary dynamics. 
We note that during inflation W = W ^ = 0 for tribrid inflation, which e.g. has the 
nice feature that it allows to combine low scale supersymmetry breaking with high 
scale inflation, as discussed in [24, 25]. 

The superfield </> contains the inflaton as scalar component. The second term in 
W generates a supersymmetric mass term for </> at the end of inflation when (H) ^ 0. 
During inflation, the inflaton potential would be flat at tree level and in global SUSY 
(with the vacuum energy during inflation given by Vo = k 2 M 4 ). The desired slope of 
the potential is generated by loop corrections (with SUSY broken during inflation due 
to the non-zero F x ) or from higher-dimensional effective operators in K. Consistency 
of the amplitude A s of the scalar perturbations and of the value of the spectral index 
n s with observations can be achieved, e.g., with the parameter choices M ~ 5 • 10 15 
GeV, k ~ 0.1, and A ~ 0.1 [26]. These values are chosen to match the set of viable 
parameter values given above for conventional SUSY hybrid inflation models. When 
we estimate the gravitino production in both classes of inflation models, we will use 
these values as benchmark point. 

After inflation, we have to consider the dynamics of three fields X, F[ and q A 
X and F[ mix almost maximally (when the SUSY breking effects after inflation are 
included) and the relevant mass eigenstates are 

x± = Vx ± H) , (2.4) 

with nearly degenerate masses m 2 ± ~ 4 k 2 M 2 . The minima are at (A") = 0, and 
(H) = M. 

The inflaton field </> has its potential minimum at (0) = 0, and this is true even 
after including SUSY breaking effects. At the minimum the inflaton mass is given 
by m 2 ~ 4A 2 M 4 /M|,, with M P being the reduced Planck mass M P = 2.4 ■ 10 18 GeV. 
As we see, it is suppressed with respect to the mass m 2 ± of the mixed A", H fields by 
M 2 /M P . Furthermore, due to the second term in W of Eq. (2.3), the \± fields can 
decay into a pair of 0 components, as we will discuss in more detail below. 

t A more general form reads W = kX (H e — M 2 ) + A H m </>", with £,m,n> 2 [22, 23]. 
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In summary, we can say that the dynamics at the end of inflation for the case of 
“conventional” SUSY hybrid inflation is governed by the oscillations and the decay 
products of the tp± field. For the case of tribrid inflation, both x± and (j) will oscillate, 
and x± will decay first (dominantly into 0 components plus some gravitinos), followed 
by the decay of the inflaton field <fi. The goal of this paper is to calculate the gravitino 
production from the decays of these fields. 

SUSY breaking sector: As we are going to discuss in the next section, the 
decay of the inflaton into gravitinos is caused by the mixing between the inflaton and 
the superfield Z which after inflation breaks SUSY via its F-term. We will assume 
for simplicity that there is only one such field, that R parity is conserved, and that 
the gravitinos have a comparatively large mass, > 30 TeV, such that they decay 
before BBN. 

Under these conditions, the gravitino decays will not spoil BBN. On the other 
hand, each gravitino will in general produce one LSP which contributes to the dark 
matter abundance. This leads to a constraint from not overproducing dark matter, or 
one may even explain the dark matter abundance from this non-thermal production. 
Since we are considering small held inflation models in this work, we will not specify 
a particular form of the Kahler potential, but rather allow for a general expansion 
in terms of helds/Mpp We will include effective operators in the Kahler potential, 
such as 7 x|A "| 4 mentioned above, or the ones which will be mentioned below. 

Let us discuss the role of the Z held after inflation: Its fermionic component 
gets eaten up by the gravitino in the super-Higgs mechanism to give the gravitino 
its mass. The scalar component Z, on the other hand, remains in the theory and 
in general also adds to the dynamics. When it is displaced from its true minimum 
at the end of inflation and if it has mass and decay rate of sizes similar to the ones 
of the gravitino, it can also lead to a severe “modulus problem” [27, 28]. It may 
dominate the universe at late times and when it decays too late it can again spoil 
the predictions of BBN or, even when it decays before BBN, its decay into two 
gravitinos would reintroduce the gravitino problem. 

However, the problem can be overcome easily with terms — ^zzlZ^/Mj, and 
— lzx\X\ ~\Z\ 2 /Mp present in the Kahler potential. 7 zz increases the mass of the 
Z-modulus held, such that m z m 3 / 2 . Furthermore, with 7 Z x > 10 the Z-modulus 
adiabatically tracks the minimum of the potential, such that only negligible amount 
of energy is transferred into Z and the modulus problem is solved [29, 30]. In the 
following we will assume suitable Kahler potential terms ^x, 7 zz and 7 Z x, such that 
we can ignore the dynamics of Z and X and that after inflation >)%> m 3 /2 
is satished. 

3 Non-thermal gravitino production in SUGRA 

As mentioned above, the gravitino can be produced either thermally or non-thermally. 
The former production takes place in the thermal bath after reheating at the end of 
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inflation, and its abundance increases with the reheat temperature T R . In addition, 
gravitinos can also be produced non-thermally, e.g. from the decay of the inflaton 
field. 

The gravitino abundance might pose a threat to inflation models: For instance, 
when the gravitinos are unstable and decay too late, they may spoil the predictions 
of BBN. However even when they decay before BBN, which is typically the case when 
they are heavier than about 30 TeV, they impose a serious constraint. Assuming R 
parity is conserved, each gravitino decay produces one LSP, which adds to the DM 
abundance, and it must not exceed D_dm^ 2 = 0.12 [2], This results in the constraint 


'mLSpYz/2 = rriLSP 


(=f) 


Plsp 

s 


< 

r^j 


Pdm 


n DM — = 4.2- 10“ 10 GeV, 
s 


(3.1) 


where Y£, 2 is the gravitino to entropy ratio produced non-thermally from the decay 
of some scalar field ip and rnisp is the mass of LSP. We have used the present critical 
density p c = 1.88 • 10” 29 h 2 g cm~ 3 , and the entropy density s = 2969 cm' 3 . For our 
discussion, we neglect the thermally produced gravitino abundance (i.e. Yy^ rmal ), 
which is justified for sufficiently small T R (i.e. typically T R -C 10 10 GeV). We will 
consider the case that DM mainly consists of thermally produced LSPs. To derive a 
bound on the gravitino abundance, we nevertheless require that the non-thermally 
produced dark matter density does not exceed Ddm^ 2 = 0.12. Furthermore we will 
neglect annihilations of the non-thermally produced LSPs, which is justified for an 
LSP which can yield the correct thermal relic density. 

Let us now turn to the non-thermal production: The part of the Lagrangian 
density which governs the decay of any scalar field ip into a pair of gravitinos is (see 
e.g. [31]) 


e~ l C = -\ e ^{G v d p ip+G z d p Z-n.c)% lv ^ a -\e G ' 2 {G^G z Z^.c )^. 
o o 

(3.2) 

Here Z is the field which breaks SUSY by its F-term (and we assume only one field 
here for simplicity), -0^ is the gravitino field and G = K + ln|VF| 2 with K and 
W being the Kahler potential and superpotential, respectively. In G v and Gz, the 
subscript denotes the derivative with respect to the field <p or Z. 

The field Z can mix with ip, and this mixing induces the decay of <p into grav¬ 
itinos, as discussed e.g. in [32, 33]. The real and the imaginary components of the 
scalar field have the same decay rates to two gravitinos, given by 



r(v> 2^3/2) ^ 


\Gf\ 2 rn% 

288vr m 2 /2 M 2 ; ’ 


(3.3) 


where |G| fl | 2 := |G0,| 2 + |e*G^| 2 (following the notation of e.g. [33]), and e* is related 
to the mixing angle between the relevant scalar field and the SUSY breaking field. 
We have assumed that the inflaton mass (after inflation) satisfies m v 7713 / 2 - In the 
above expression, the only model dependent parameter is and for the considered 
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case, rriz m 3 / 2 , it is given (in Planck units) by 


\Gf r ^ 


r- m 7 
v39u>z 5 " 


m- 


+ 


V3 (V,Gz)^ 

mt 


+ 


3g<pzz 


m 3 / 2 


m„ 


(3.4) 


Here g t j denotes derivatives of the Kahler potential with respect to the field (pi and 
<fj. The first term in Eq. (3.4) thus contributes whenever there is non-zero mixing 
between ip and Z in the Kahler metric. When ip is charged differently then Z (or Z) 
under some symmetry (including R symmetry), then this term is always proportional 
to (<p). Regarding the second term, we note that | ~ (<p) when the Kahler 

potential contains a canonical term for (p , K = \<p\ 2 + ... [33]. The third term 

proportional to g^zz can contribute significantly to the decay rate as long as the 
SUSY breaking field is uncharged such that terms like a\ip\-(Z 2 + Z )/Mp exist in 
K. The contribution to the decay rate from this term is then also proportional to 
the field vev (ip). 

In the following, we will apply these results to the fields which are relevant 
after inflation, which in “conventional” SUSY hybrid and tribrid inflation models 

are (primarily) the inflaton field and the waterfall field. Such fields can acquire non¬ 

zero values of G| ff from the terms in Eq. (3.4), and then produce gravitinos when 
they decay. In general all three of these terms may contribute, however in order to 
keep the discussion simple we will focus on the effect of the last term of Eq. (3.4) and 
compare the resulting gravitino production in hybrid and tribrid inflation models. 
We choose this case as an example since when the third term is non-zero, it typically 
provides the dominant contribution to G| ff . 


4 Non-thermal gravitino production in hybrid inflation 


As mentioned above, we will focus on the case that the contribution from the third 
term in Eq. (3.4) dominates the decay rates to gravitinos. More explicitly, we consider 
as example the effects of a Kahler potential term 5K = a\H\ 2 (Z 2 + Z )/Mp , which 
exists when Z is uncharged. We then obtain 

(H) my 2 


|Gf + |^3a- 


M P m^ ± ’ 

and the decay rates of the mass eigenstates to gravitinos are given by* 


* 9q 2 (g) 2 m; 

V 2 288ir M% 


3 

<t>± 


(4.1) 


(4.2) 


The (p± fields decay mainly to the SM particles and their superpartners (with a 
total decay rate Tf ot ), giving rise to the reheat temperature 

-1/4 


n = 


r LMp 


n 2 g* 

10 


(4.3) 


*We note that TgU = Tg, 2 =: I'W. We will nevertheless continue to use m^ ± to distinguish it 
from the post-inflationary inflaton mass in tribrid inflation. 
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where g* is the effective number of relativistic degrees of freedom at the temperature 
T r . Counting supersymmetric degrees of freedom, we take g* ~ 228. Neglecting 
gravitino production in the thermal bath, the (non-thermally produced) gravitino 
number density from 0 decays K^ 2 , normalized to entropy density, is given by 



r 


~ 2 


V O rp(p 

3/2 3 J-R 


r L 4% 


3 / 7T 2 g* \ 2 r 3/2 Mp 

2 V 10 / ' 


(4.4) 


Using Eq. (4.2), this leads to 


m LSP Yy 2 ~ kV f 


<#> 


10 16 GeV 


Af 


5 • 10 15 GeV 


m LS p \ 

100 GeV/ 


/10 7 GeV\ 

V / 


10 GeV, 
(4.5) 


where we have used m^ ± ~ 4 r 2 M 2 . Assuming a ~ C>(1), and for typical values 
of (H) ~ 10 16 GeV, M ~ 5 • 10 15 GeV, k ~ 0.1 [19], we see immediately that 
the estimate violates the observational constraint of Eq. (3.1) for a wide range of 
parameters. Unless the higher order terms in the Kahler potential are highly sup¬ 
pressed, say a ~ 10 -5 , the gravitino production is problematic for “conventional” 
SUSY hybrid inflation models “. While the thermal production of gravitinos can be 
reduced by lowering Tf 0 this would increase the non-thermal gravitino production in 
“conventional” SUSY hybrid inflation. 


5 Non-thermal gravitino production in tribrid inflation 

In this section we estimate how many gravitinos are produced non-thermally in tribrid 
inflation, i.e. from the decay products of the inflaton 0, the waterfall field H , and the 
“driving field” X. As it has been mentioned earlier, at the end of inflation H and 
X mix maximally to form the mass eigenstates x± (with nearly degenerate masses). 
Therefore, the dynamics of x± an d 0 is relevant for calculating the non-thermally 
produced gravitino abundance. 


5.1 Decay rates in gravitinos 

As it has been noted in Eq. (3.3), the decay rate to gravitinos from a scalar field 
depends on (given by Eq. (3.4). For tribrid inflation G{f vanishes since (0) = 0 
in the minimum after inflation, and only is nonzero. We will again assume 
that the third term in Eq. (3.4)) dominates and focus on the effects from a term 
SK = a\H\ 2 (Z 2 + Z~)/Mp[ in the Kahler potential. In total, we obtain 


| neff 

\<? 


0, 


eff I 

X± I 


3a 


ix) "U /2 


Mp m 


x± 


(5.1) 

(5.2) 


§One way to avoid the gravitino problem in “conventional” SUSY hybrid inflation arising from 
the third term of Eq. (3.4) is by charging Z under a symmetry such that g$ zz vanishes. Such a 
case has been studied, e.g., in [34]. 
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This means that only the decays of x± will produce gravitinos, while the gravitino 
production from inflaton decays will be practically zero. 


5.2 Estimation of the produced gravitino number density 


To estimate the gravitino number density we will make various simplifying approxi¬ 
mations: We will begin our considerations when the fields x± decay. We will assume 
that at this time the amplitude of the 0 field oscillations is sufficiently small and that 
the energy density in the coherent oscillations simply redshifts like non-relativistic 
matter. In addition, we will divide the post-inflationary epoch in either radiation 
dominated or matter dominated phases. These approximations will allow us to de¬ 
rive rather simple analytical formulae. Of course, to calculate gravitino production 
in an explicit model numerically, one would drop these assumptions. 

We will furthermore assume that T} ot > rf ot , i.e. that x± decays earlier than 
the inflaton field 0, and also that m x ± > 2 m^. This is generically true for most pa¬ 
rameter choices in tribrid inflation consistent with observations (in particular for the 
parameter sets we use here as examples). With m x± > 2m^, x± decays dominantly 
in two inflatons or its SUSY partners with total decay width approximately given by 
(see e.g. [10] 


1 tot — 


A 2 df 3 
8irMp ’ 


(5.3) 


One important parameter for our analytic considerations will be the ratio of the 
energy densities stored in the x± field (in terms of homogenous scalar field oscilla¬ 
tions) and 0 field when x± decays, i.e. 


« = Px/P<fi > ( 5 - 4 ) 

where the tilde denotes the time of x± decay. In the following we will treat cc as a 
free parameter, as it depends on the details of the model and as its calculation would 
require the simulation of the preheating stage after tribrid inflation, which is beyond 
the scope of this paper. 

At the time of x± decay all the decay products are relativistic and their energy 
densities therefore redshift as radiation. We can write for the energy densities p x = 
prad _|_ p ^ w ] iere p 3 ^ 2 w pp prad ( ] eno ti n g a ll other radiation components 

from x± decay. 

We will now estimate the gravitino production for different regimes of a: 

• If a < 1, the radiation energy density produced from x± decay is smaller than 
the energy density of the coherently oscillating 0 field, which means that in 
our approximation we will treat the universe as matter dominated. It then 
remains matter dominated untill the 0 field decays into radiation, and at this 
time we calculate the reheat temperature. During the phase where the 0 field 
dominates, the gravitinos from x± decay, which are treated as radiation, get 
diluted. 

Un addition, \± decays (with lower rates) also to gravitinos as well as to MSSM particles. 
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• On the other hand, if a > y r x ot /Tf ot (as we will discuss below), the universe 
is (mainly) radiation dominated already after the time of x± decay (which is 
the case for a > 1), and it will remain radiation dominated until the time of (p 
decay. The reheat temperature is then calculated at the time of x± decay. 


• For intermediate values of a , the universe is radiation dominated directly after 
X± decay, but becomes matter dominated before (p decay. In this case, the 
reheat temperature is calculated at the time of (p decay. 


Case ex. < 1: As mentioned above, in our approximation we will treat the universe 
as matter dominated in the epoch between x± decay and (p decay. The decay of the 
inflaton field <p dilutes the gravitino abundance produced from x± decay, and from 
<p decays no additional gravitinos are produced. The resulting gravitino abundance 
can be estimated as follows: 

The number density of gravitinos at the time of x± decay is given by 

T x 

^ 3/2 = fi x ^r ■ 2 , ( 5 - 5 ) 

where the factor 2 comes from the fact that x± decays into two gravitinos. For a < 1, 
the Hubble constant at the time of <p decay, FT, can be related to the Hubble constant 
at the time of x± decay, FF, by 


FF 

FF 




So we obtain for the number density at the time of <p decay: 


™3/2 



(5.6) 


(5.7) 


After <p decay, the gravitino number density normalised to the entropy density, Y :i / 2 
is thus given by 


Y 3/2 = 


™3/2 



(5.8) 


where p x = m x n x has been used. Assuming thermal equilibrium at the end of (p 
decay for simplicity, we use 


£ trad, %rad o 

Ptot Ptf> P x Prad ^ rp 

~~ — ——I-— — — T J-R, 

s s s s 4 


(5.9) 


which leads to 

? = i Ts (i + “ ( r «/r?oi) _2/3 ) ■ (s.io) 

Tf t is the reheat temperature that has contributions from both <p and x decays. 
Plugging in this relation gives us the following expression for the gravitino yield: 


* 3/2 


2 afrjA 

m x V T tot ) 


3 T r 
4 


^1 + OL 



(5.11) 
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The second term in the bracket becomes negligible when either a <C 1, i.e. when 
most of the energy density after inflation is carried by the inflaton field 0 (in which 
case there is an overall suppression of Y 3/2 by a factor a), or when Tf ot Tf ot , i.e. 
when 0 decays much later than x±- Under one of these conditions, the above formula 
reduces to the simple form 


Interestingly, and in contrast to the case of “conventional” hybrid inflation, the 
gravitino yield is proportional to the reheat temperature Tr. Therefore by lowering 
Tr, which corresponds to the case that the inflaton field 0 decays comparatively 
late, we can strongly suppress the non-thermally produced gravitino density V 3 / 2 in 
tribrid inflation. We will now plug in values for the model parameters (as we did for 
“conventional” SUSY hybrid inflation) to check that for tribrid inflation the gravitino 
production can be sufficiently suppressed to avoid the “gravitino problem”. 

For a rough estimate of the gravitino yield, we use T* ot — T^, neglecting the 
contribution from the decays to gravitinos and to Standard Model (MSSM) particles. 
The total decay rate of x is then given by Eq. (5.3). Note that this is a conservative 
approximation, and will provide an estimate for an upper bound on Y 3 / 2 . Using 
Eq. (5.2), the gravitino yield is given by 


^ _ 3 cm 2 mKx) 2 ^ 

3/2 8K\ 2 M 3 My R 


(5.13) 


Using the post-inflationary mass for the x± fields, m 2 ± ~ An 2 AT 2 , as in Eq. (4.5) 
for the case of “conventional” SUSY hybrid inflation, we find the corresponding 
expression for tribrid inflation 


rriLSpY ; 3 /2 - 


Aaa 2 n 


(x) 


X 2 \ 10 16 GeV 


5-10 15 GeV\/ m LS p \ 


M 


\ ( m LSP \ ( T r A 12 
) VlOOGeVAlO 7 GeV/ 


(5.14) 

Plugging in the typical values for tribrid inflation (x) ~ 10 16 GeV, M — 5 • 10 15 GeV, 
hi = 0.1, A = 0.1 [ 6 , 26], we see that the abundance is well bellow the observational 
limit of Eq. (3.1) for m R sp ~ 100 GeV and for the example reheat temperature 
T r < 10 7 GeV. 


Note that for reheat temperatures around 10 7 GeV the thermal production of 
gravitinos can indeed be neglected. We like to emphasize that, in contrast to the 
non-thermal gravitino production in standard SUSY hybrid inflation, the gravitino 
abundance in tribrid inflation is proportional to the reheat temperature T R . Thus by 
lowering Tr, in tribrid inflation models one is also reducing the non-thermal gravitino 
production. 


Case 1 < a < \/r* ,/lt, : 

If 1 < a < the decay of the \ field would produce a radiation 

dominated era which would be followed by a matter dominated phase with coherent 
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Figure 1. The diagram illustrates the different times and epochs used in the calculation 
of our estimates for the gravitino number density Y 3/2 in tribrid inflation. For a —> 1, the 
time t*, where radiation and matter energy densities are equal (shown by the green line 
in the diagram) is moved further towards the left and merges with t, the time when x± 

decays. For a —>■ \Jrf 0f /rf 0 f, t* moves towards the right and merges with the time t when 
0 decays. 


oscillations of the 0 field. The time at which radiation and matter densities become 
equal, will be denoted by (cf. Fig. 1 for a schematic illustration). Using H oc a -2 
for the radiation dominated epoch and H oc a -3 / 2 for the time of matter domination, 
we arrive at 

a,/a = « 4 / 3 (rL/rf ot ) , (5.15) 

where H * = r^(a/a *) 2 has been used. This translates to the following condition for 
having both, a radiation dominated and a matter dominated phase between x± an d 
0 decay: 

1 < a < 0n„,/r? oi . (5.16) 

Again, tracking the energy densities during different epochs we arrive at 


a/a = a 




(5.17) 


where a = p x /p<f> = a*/a has been used. Following similar steps as for the case 
a < 1 , we find that the gravitino number density Y 3/2 can be estimated as 


>3/2 


2a AV\ 

m x \ r L ) 


370 

4 




(5.18) 


We see that (as it should be) for a = 1 the above formula matches to the one in 
(5.11). 


Case a > 

Ill the limit a —1 (r(,i/n(. the time t* where radiation and matter energy den- 
sities are equal merges with the time t when 0 decays. This means that there will be 
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no phase of intermediate matter domination and thus in our approximation we treat 
the universe as radiation dominated from the time of x± decay on. Consequently, 
in our approximation, there will be no further dilution of the gravitinos produced 
from x± decay. Furthermore, the reheat temperature is now defined at the time of 
X± decay. We therefore expect that in some approximation one can calculate the 
gravitino abundance as in “conventional” SUSY hybrid inflation from x± decay, i.e. 
using the formula from Eq. (4.4) with 0 replaced by y: 


Y 3 /2 


rL 4 m x± 


(5.19) 


To see that is indeed true, let us consider the limiting case a = \jY} ot /Tf ot in more 
detail: For this value, = p r x d , and the reheat temperature, T R = | (p™ d + p™ d )/s 
can be written as T R = |(2 p™ d )/s =: 2T^, defining : = \p™ d /s as the reheat 

temperature calculated from 0 decay alone (ignoring p r x d ). Defining analogously 
T r '■= \p r x d /s (at the time of y± decay), using 

Tr/t% = \frijrt , (5.20) 


and plugging in \JYf ot /Tf ot = a, we indeed verify that Eq. (5.18) reduces to Eq. (5.19). 
Note that using T* { in Eq. (5.19) means that the contribution p^ is neglected when 
calculating the reheat temperature at the time of x± decay (which is the approxima¬ 
tion in which both formulae match). 

Let us note at this point that F f Y ot /rf oi is actually a quite large number. Con¬ 
sidering again the limiting case \jYf ot /Tf ot = a where T ^ = \T r , and plugging in 
the above used example values for the tribrid inflation model parameters (y) ~ 10 16 
GeV, M = 5 • 10 15 GeV, k = 0.1, A = 0.1 [6, 26], we obtain (with Tf ot from Eq. (5.3) 
and r f ot from Eq. (4.3)): 


2M 3 \ 2 k (n yyj in 9 ( io 7 GeV V 

r t t rM p svt V io ) V t r ) 


(5.21) 


This implies in particular that a is very large, a = \JY} ot /Tf ot ~ 5 ■ 10 4 for T R ~ 10' 
GeV, which means that it is a very good approximation to neglect p$ when calculating 
the reheat temperature at the time of x± decay (since p x = ap^ p ^). On the other 

hand, we note that we do not expect that such a large a, where a > \jYf ot /Tf ot , can 
be realised. Therefore, the non-thermally produced gravitino abundance in tribrid 
inflation should be generically suppressed. 


6 Conclusions 

In this paper, we have investigated non-thermal gravitino production after tribrid 
inflation. Tribrid inflation is a variant of hybrid inflation where the inflaton resides in 
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the matter sector of the theory (cf. Sec. 2). We find that in contrast to conventional 
supersymmetric hybrid inflation, where non-thermal gravitino production imposes 
severe constraints on the inflationary model, the “non-thermal gravitino problem” is 
generic-ally absent in models of tribrid inflation. 

As in “conventional” SUSY hybrid inflation, at the end of tribrid inflation the 
dynamics is governed by oscillations and decays of the inflaton field and the waterfall 
field. However, in contrast to “conventional” SUSY hybrid inflation, these fields are 
not degenerate in mass but (in the considered tribrid inflation setup) the inflaton 
field is much lighter than the waterfall field. While the decays of the waterfall field in 
tribrid inflation (and of both fields in “conventional” SUSY hybrid inflation) produce 
gravitinos as discussed in Sections 4 and 5, the inflaton field in tribrid inflation has 
a negligible decay rate into gravitinos. The two main effects which allow to easily 
solve the “non-thermal gravitino problem” in tribrid inflation are: 

• With the inflaton being lighter than the waterfall field, in tribrid inflation the 
latter automatically has a second decay channel with a much larger rate than 
for the decay into gravitinos. This leads to a strongly suppressed gravitino 
production from the decays of the waterfall field. 

• The inflaton in the considered tribrid inflation setup generically decays later 
than the waterfall field. While it does not produce gravitinos from its decays 
(cf. subsection 5.1), it even dilutes the gravitino population produced earlier 
from the waterfall field decays. The later the inflaton decays, the smaller the 
reheat temperature Tr, and the stronger is the dilution effect. 

To quantify how much the “non-thermal gravitino problem” is relaxed, in sub¬ 
section 5.2 we calculated estimates for different regimes of a, which parameterizes 
the ratio of the energy densities in the waterfall field and in the inflaton field after 
inflation. Comparing the produced gravitino abundances with conventional F-term 
SUSY hybrid inflation, we find that the overall amount of produced gravitinos is in 
general much lower (cf. e.g. Eqs. (4.5) vs. (5.14)), and also that there is a different 
dependence on the model parameters. For instance, in tribrid inflation a lower Tr al¬ 
lows to suppress the produced gravitino abundance, whereas in conventional F-term 
SUSY hybrid inflation, the gravitino yield is proportional to 1 /Tr. 

In summary, models of tribrid inflation not only offer the attractive possibility to 
identify the inflaton with a field (or a H-flat direction of fields) from the matter sector 
of the theory. As we showed in this letter, they also allow to solve the “non-thermal 
gravitino problem” easily due to some built-in model properties. 
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